Background: 


« Astronauts may lose up to 9% of load- 
bearing bone density per month in 
Spaceflight! 

Lower chance of fracture in space due 
to lower loads (OG) 2 


Higher loads on Earth (1G) result in a higher potential 
for fracture due to lowered bone density when 
astronauts return to Earth 


Computational bone strength model can be used to 
assess bone fracture risk for astronauts 


Objective: 

Develop and test an open-source computational bone 
strength model for acceptable performance in the 
assessment of pre-flight and post-flight astronaut bone 
Strength studies. 


Open Source Advantage : 

¢ Publicly published with a community collaborative 
mindset, where others are encouraged to view and 
contribute to the code to advance development 

¢ Allows for expanded future development and input 
from a large community of experts 


Hypothesis: 


¢ Combine existing open-source software with our own 
scripts (Python) 
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Performing image segmentation through 
Python 
Image Segmentation: 
« Python script imports CT scans with visualization toolkit 
(VTK) 
¢ Library allows for import of many popular medical 
image formats 
¢ Script translates pixel values to Hounstfield values 
using metadata in original CT scans 
¢ Script isolates bone from medical images with 
thresholding based on Hounsfield values 


¢ Final image is binary representation of bone regions 
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3D Model Construction: 
on 


¢« Python's VTK toolkit includes a 
Marching Cubes algorithm python 


¢ Creates a 2D surface mesh 
from binary segmented 
regions? 

¢ Also smooths mesh and 
removes unnecessary 
triangles 


¢ Blender* used to repair mesh and 
isolate any areas of interest 

¢ 2D surface mesh was recreated into 
a 3D volume mesh with Gmsh» 
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Flowchart showing 3D 
model creation process 
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Isolated vertebrae disc from processed skull mesh 


Material Properties Assignment: 
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Rescale FEBio 
mesh to CT 
mesh 
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| Current process of assigning material properties to mesh 
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¢ Original CT scans’ Hounsfield 
values extrapolated into 
densities and Youngs modulus 
using Keneko et al.'s © prior bone 
ash testing 
¢ Translates Hounstield 
value to bone ash density 
¢ Extrapolates Youngs 
modulus from bone ash 
density 
¢ Python script writes material 
properties to an FEBio XML file 


Skull model with dynamically 7 
for easy import 


applied material properties 
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OCT Dengity, p57 (myer) 
Figures from Keneko et al. ¢ showing relationship between Hounsfield value 
ash density, and Youngs modulus 
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In-House Developed Interface: 


B.| NASA Glenn DICOMesher ? x 


¢ A Graphical User 
Interface (GUI) combines 
process of image 
segmentation and 2D 
mesh creation into a 
centralized tool 


Load DICOM Stack Image Segmentation 2D Mesh Creation 


Select Directory of DICOM Images 


Browse ColUsessAbbishelthel Google Drive NASADO1SUmaging 0S 


ee # of Images 
coat 361 


View Images 


[ Viewimage sd 
“™ —e User can select stack of 
m python medical images to import 


and can view stack in 
three dimensions 


¢ Users can perform image 
segmentation using their 
inputted threshold value 

¢ Exports a 2D mesh for 
next processing step In 
Blender and Gmsh 


Graphical User Interface for tool 


Contour map of effective stress (left) and pressure (right) on vertebrae model 


Finite Element Analysis (FEA): 


¢ Finite Element Analysis performed through FEBio’ suite 


¢ Software allows for the graphical fixing of points, 
defining of loads and boundary conditions 


¢ Allows for graphical viewing of end results 


Conclusions: 
¢ No straightforward method to implement existing 
open-source software into desired product 
¢ Acombination of various open source software along 
with self-developed scripts was needed to complete 
the segmentation, 3D construction, and FEA analysis 
tasks 


Future Work: 
¢ Need to design and run a selection of test cases to 
validate our method, including a full end-to-end 
simulation 
¢ Extend further aspects of tool into interface, allowing 
for Tull integration of method into a single location 
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